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. Abstract. An extensive program for the calculation of galactic cosmic-ray prop- 

P\J ' agation has been developed. This is a continuation of the work described in [1] . 

The main motivation for developing this code [2] is the prediction of diffuse Galac- 
tic gamma rays for comparison with data from the CGRO instruments egret, 
COMPTEL, and OSSE. The basic spatial propagation mechanisms are (momentum- 
dependent) diffusion, convection, while in momentum space energy loss and dif- 
^SJ ■ fusive reacceleration are treated. Primary and secondary nucleons, primary and 

On ' secondary electrons, and secondary positrons are included. Fragmentation and 

. energy losses are computed using realistic distributions for the interstellar gas 

' and radiation fields. 

o 

c/3 . We are developing a model which aims to reproduce self-consistent ly ob- 

^ ' servational data of many kinds related to cosmic-ray origin and propagation: 

direct measurements of nuclei, electrons and positrons, gamma rays, and syn- 
^ ■ chrotron radiation. These data provide many independent constraints on any 

model and our approach is able to take advantage of this since it must be 
consistent with all types of observation. We emphasize also the use of real- 
istic astrophysical input (e.g., for the gas distribution) as well as theoretical 
developments (e.g., reacceleration). The code is sufficiently general that new 
physical effects can be introduced as required. The basic procedure is first to 
obtain a set of propagation parameters which reproduce the cosmic ray B/C 
ratio, and the spectrum of secondary positrons; the same propagation con- 
ditions are then applied to primary electrons. Gamma-ray and synchrotron 
emission are then evaluated. Models both with and without reacceleration are 
considered. The models are three dimensional with cylindrical symmetry in 
the Galaxy, the basic coordinates being {R,z,p) where R is Galactocentric 
radius, z is the distance from the Galactic plane and p is the total particle 
momentum. The numerical solution of the transport equation is based on a 
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Crank-Nicholson implicit second-order scheme. In the models the propagation 
region is bounded \)Y z — Zh beyond which free escape is assumed. A value 
= 3 kpc has been adopted since this is within the range which is consis- 
tent with studies of ^°Be/Be and synchrotron radiation. For a given Zh the 
diffusion coefficient as a function of momentum is determined by B/C for the 
case of no reacceleration; with reacceleration on the other hand it is the reac- 
celeration strength (related to the Alfven speed va) which is determined by 
B/C. Reacceleration provides a natural mechanism to reproduce the B/C ratio 
without an ad- hoc form for the diffusion coefficient [3,4] . The spatial diffusion 
coefficient for the case without reacceleration is D = PDq below rigidity po, 
(3Dq{p/ pqY above rigidity po, where (5 is the particle speed. The spatial diffu- 
sion coefficient with reacceleration assumes a Kolmogorov spectrum of weak 
MHD turbulence so = j3Do(p/ poY with 6=1/3 for all rigidities. For this 
case the momentum-space diffusion coefficient is related to the spatial one 
[4]. The injection spectrum of nucleons is assumed to be a power law in mo- 
mentum. The interstellar hydrogen distribution uses HI and CO surveys and 
information on the ionized component; the flelium fraction of the gas is taken 
as 0.11 by number. The interstellar radiation field for inverse Compton losses is 
based on stellar population models and IRAS and COBE data, plus the cosmic 
microwave background. Energy losses for electrons by ionization. Coulomb, 
bremsstrahlung, inverse Compton and synchrotron are included, and for nu- 
cleons by ionization and Coulomb interactions. The distribution of cosmic-ray 
sources is chosen to reproduce the cosmic-ray distribution determined by anal- 
ysis of EGRET gamma-ray data [5] . The bremsstrahlung and inverse Compton 
gamma rays are computed self-consistently from the gas and radiation fields 
used for the propagation. The 7r°-decay gamma rays are calculated explicitly 
from the proton and Helium spectra using [6]. The secondary nucleon and 
secondary source functions are computed from the propagated primary 
distribution and the gas distribution, and the anisotropic distributions of 
in the //^ system was taken into account [7]. 



ILLUSTRATIVE RESULTS 

Some results obtained are shown in the Figures. The energy dependence 
of the B/C ratio, and local proton and Helium spectra are shown in Figs. 1 
and 2. The spectrum of primary electrons is shown in Fig. 6. The adopted 
electron injection spectrum has a power law index -2.1 up to 10 GeV; this 
is chosen using the constraints from synchrotron and from gamma rays. The 
electron spectrum is consistent with the direct measurements around 10 GeV 
where solar modulation is small and it also satisfies the constraints from J'^ , . 

e + e+ 

Above 10 GeV a break is required in the injection spectrum to at least -2.4 for 
agreement with direct measurements (which may however not be necessary if 
local sources dominate the directly measured high-energy electron spectrum). 
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The synchrotron spectrum at high Galactic latitudes (Fig. 3) is important 
since its shape constrains the shape of the 1-10 GeV electron spectrum. An 

injection index -2.1 (without reacceleration) is the steepest which is allowed 
by the radio data over the range 38 to 1420 MHz. As illustrated, an index -2 A 
as often used (e.g., [8]) gives a synchrotron spectrum which is too steep. 

The modelled gamma-ray spectrum for the inner Galaxy, shown here for the 
case of no reacceleration (Fig. 4), fits well the comptel [9] and EGRET [5] 
data between 1 MeV and 1 GeV beyond which there is the well-known excess 
not accounted for by 7r°-decay with the standard nucleon spectrum [10,11]. 
The electron spectrum used here is increased by a factor 2 over that shown 
in Fig. 6; this is required to reproduce well the observed gamma intensities. 
It is still within the range allowed by the positron fraction (see below), and 
the shape accords well with the flatter electron injection spectrum required 
by synchrotron data (Fig. 3). Inverse Compton is dominant below 10 MeV, 
bremsstrahlung becomes important for 3-200 MeV. The lower bremsstrahlung 
combined with 7r°-decay leads to a good fit to the flat spectrum observed in 
this range, in contrast to previous attempts to model the spectrum with a 
steeper bremsstrahlung spectrum [8]. The 10-30 MeV 7-ray longitude profile 
(Fig. 5) at low latitudes from this model can be compared with that from 
COMPTEL [9]. The e+ spectrum (Fig. 6) and the positron fraction (Fig. 7) 
agree well with the most recent data compilation for 0.1-10 GeV [7]. 

This study indicates that it is possible to construct a model satisfying a 
wide range of observational constraints and provides a basis for future devel- 
opments. 

More details can be found on http://www.gamma.rnpe-garching.rnpg. de/'^aws/aws.htnil 
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FIGURE 1. The energy dependence 
of the B/C ratio can be reproduced 
with Do = 2.0x10^^ cxn^s'^, 6 = 0.6, 
po = 3 GV/c without reaccelera- 
tion (thick hue) and Dq = 4.2 xlO^^ 
cm^s^^j VA = 20 km s^^ with reac- 
celeration (thin hue). Dashed hnes are 
modulated to 500 MV. Data from [12]. 




kinetic erte-rgy, MeV/nticleo-n kinetic energy /nucleon, MeV 



FIGURE 2. Left panel: the local proton spectrum for injection index 2.15 (thin solid 
line), 2.25 (thick solid line) without and with reacceleration respectively, compared with the 
measured 'interstellar' spectrum (dashed [13] and dashed-dot [11] lines). Right panel: the 
Helium spectrum with injection index 2.25 (thin solid line), 2.45 (thick solid line) without 
and with reacceleration respectively, compared with the measured 'interstellar' spectrum 
(dashed [13] and dashed-dot [14] lines). 




FIGURE 3. The synchrotron spectrum at FIGURE 4. The 7-ray spectrum for the 

the NGP, and predictions for electron injection inner Galaxy, 330° < I < 30°, |6| < 5°. 
indices -2.1 (dashed line) and -2.4 (solid line), egret data [5], comptel data [9]. 
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FIGURE 5. Longitude distribution of gamma rays in energy range 10-30 MeV. Histogram: 
COMPTEL [9], dashed: inverse Compton, thin line: bremsstrahlung, thick hne: total. 
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FIGURE 6. Spectra of secondary e^, and of primary electrons. Full lines: our model with 
no reacceleration (left) and with reacceleration (right). Dash-dotted lines by Protheroe [15]: 
lower is his leaky-box prediction for e"*", upper is his adopted electron spectrum. 
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FIGURE 7. Positron fraction for model with no reacceleration (left panel) and with reac- 
celeration (right panel). Dashed and dash-dotted lines by Protheroe [15]: his predictions of 
the leaky-box and diffusive halo models respectively. The data collection is taken from [16]. 



